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Stellar models for very low mass main sequence stars: the 
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ABSTRACT 

We present Very Low Mass stellar models as computed including non-grey model atmospheres for selected assumptions 
about the star metallicities. The role of atmospheres is discussed and the models are compared with models based on the 
Eddington approximation and with similar models appeared in the recent literature. Theoretical predictions concerning both 
the HR diagram location and the mass- luminosity relation are presented and discussed in terms of expectations in selected 
photometric bands. Comparison with available observational data concerning both galactic globular clusters and dwarfs in 
the solar neighborhood reveals a satisfactory agreement together with the existence of some residual mismatches. 

Key words: stars: evolution - stars: low-mass, brown dwarfs - stars: population I 
and II - globular clusters: general 



1 INTRODUCTION 

The long-standing theoretical interest in Very Low Mass 
(VLM) Main Sequence stars has been recently rejuvenated 
according to the increasing amount of VLM objects observed 
with the Hubble Space Telescope as well as thanks to CCD 
parallax determinations (Monet et al. 1992, Dahn et al. 
1995, Tinney 1996) which are increasing the amount of ab- 
solute magnitudes available for nearby dwarfs. This interest 
is further enhanced by the suggestion that an appreciable 
fraction of the baryonic mass in most galaxies could be in 
the form of VLM stars and Brown Dwarfs. To have light on 
such an observational scenario one needs reliable theoretical 
predictions about VLM stellar structures and, in particular, 
accurate mass-luminosity relations allowing the evaluation 
of reliable mass functions and, in turn, reliable estimates of 
the mass density of VLM stars in the Galaxy. 

The temperature scale of M dwarfs has been for long 
time an unsettled problem and it is still a key ingredient 
for understanding the location of VLM stars in the HR dia- 
gram (Bessel 1995, Leggett et al. 1996). Since the pioneer- 
ing works by Limber (1958), Hayashi & Nakano (1963) and 
Ezer & Cameron (1967), numerous investigators (D'Antona 
1987; VandenBerg et al. 1988; Burrows, Hubbard & Lu- 
nine 1989; Dorman, Nelson & Chau 1989; Burrows et al. 
1993; D'Antona & Mazzitelh 1994, 1996; Baraffe et al. 1995; 
Baraffe & Chabrier 1995; Baraffe & Chabrier 1996; Alexan- 
der et al. 1997; Kroupa & Tout 1997) have already de- 
voted their attention to the structural properties and the 
evolutionary behavior of VLM stars. As early recognized, 
one knows that convection in VLM stars is very efficient 
throughout all the structure and that a VLM structure is 
very nearly adiabatic. As a consequence, theoretical predic- 
tions for the interior of such peculiar objects do not depend 



on radiative opacities, nor on the choice of the mixing length 
parameter governing the superadiabatic convection. 

In spite of these advantages, VLM stellar models criti- 
cally depend on the evaluations of both opacity and equa- 
tion of state (EOS) for a low temperature, high density gas, 
where molecules, grains and non-ideal gas effects play a rel- 
evant role. Concerning radiative opacities, large efforts have 
been recently devoted to properly include in opacity compu- 
tations the effects of molecules and grains, with particular 
emphasis on the treatment of the opacity due to the H2O 
and TiO molecules (see, e.g., Dorman et al. 1989, Allard 
& Hauschildt 1994, 1995a, 1995b, Bessel 1995). Low tem- 
perature opacity tables presented by Alexander & Ferguson 
(1994) and Alexander (1994) represent a significant improve- 
ment of previous evaluations given by Cox & Tabor (1976), 
Alexander, Johnson & Rypma (1983), and Kurucz (1992). 

Substantial progresses (Saumon 1994) in the EOS for 
VLM stars are due to Saumon, Chabrier & Van Horn (1995) 
(but see also Saumon & Chabrier 1992), who provided an 
EOS for dense and cool matter based on a detailed de- 
scription of the physics at work. However, the theory of 
VLM structures requires another critical ingredient, i.e., the 
boundary conditions for the inner stellar structure, as given 
by a suitable treatment of stellar atmospheres. Grey model 
atmospheres usually assume that the optically thin region 
lying above the photosphere is not affected by convection. 
However, a similar assumption is not well-grounded for sev- 
eral VLM atmospheres, since atmospheric opacities can be 
very large due to the presence of molecules whereas the adi- 
abatic gradient is small due to the dissociation of H2. As 
a consequence the atmosphere can be unstable against con- 
vection also at very small optical depth (r). Moreover, non- 
grey effects in the stellar atmospheres are often neglected, 
though, theoretical evidences (Saumon et al. 1994) suggest 
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that such effects could play a significant role in governing 
the atmospheric structure. 

In the last decade, significant improvements in model 

atmosphere for late M dwarfs have been presented by Al- 
lard (1990), Kui (1991), Brett & Plez (1993), Saumon et 
al. (1994), AUard & Hauschildt (1995) and Brett (1995a, 
1995b). Similax models, and in particular the models pro- 
vided by Saumon et al. (1994) for zero metallicity mixtures 
and by AUard & Hauschildt (1995) and Brett (1995a,1995b) 
for finite metallicitics, represent a substantial progress of 
our knowledge in the field. Nevertheless, Bessell (1995) has 
pointed out the still existing differences between different 
models, as due to the different evaluations of the contri- 
bution to opacity of H[20 and TiO molecules and to the 
different opacity averaging technique adopted in the various 
works. 

Prom am evolutionary point of view, the need for ac- 
curate model atmospheres in computing VLM stellar mod- 
els has been firstly stressed by Burrows ct al. (1993) and 
more recently reinforced by Baraffe et al. (1995), Baraffe 
& Chabrier (1996), Chabrier, Baraffe & Plez (1996), Mera, 
Chabrier & Baraffe (1996). In previous works (Alexander 
et al. 1997, hereinafter Paper I; Brocato et al. 1997a), 
we have presented a theoretical approach to the evolution- 
ary behavior of VLM stars, showing that the adoption of 
the most updated equation of state (Saumon et al. 1995) 
and low temperature opacity (Alexander & Ferguson 1994), 
but still relying on an approximate treatment of the stel- 
lar atmosphere, allows a rather satisfactory agreement be- 
tween observation and theoretical predictions for the Color- 
Magnitude (CM) diagrams and the mass-luminosity relation 
of both metal poor and solar metallicity VLM objects. We 
thus concluded Paper I with the statement that " the use 
of a T(t) relation in computing stellar models has to be re- 
garded as a first order but not-too-bad approximation to the 
expected evolutionary behavior". In this paper we will go 
deeper in that matter, discussing VLM stellar models includ- 
ing updated outer boundary conditions for various selected 
assumptions about the star metallicities. 

The layout of this paper is as follows. Next section 
will provide some general informations about the models, 
with particular attention to the adopted grid of model at- 
mospheres. In §3, models computed adopting model atmo- 
spheres will be compared with similar models based on the 
Eddington approximation and with stellar models already 
appeared in the literature. Section 4 presents the compari- 
son between observational data and theoretical results. Con- 
clusions will follow in §5. 

2 BOUNDARY CONDITIONS. 

Models presented in this paper adopt the same EOS 
and the same opacity evaluations as in Paper I. The main 
difference with models in Paper I is the different approach 

adopted for deriving the outer boundary conditions, i.e., 
temperature and pressure of the gas at the basis of the at- 
mosphere. In Paper I the atmospheric integrations were 
performed adopting the Krishna-Swamy (1966) solar scaled 
T(r) formula until reaching the optical depth t = 2/3 or, 
alternatively, until the onset of convection, whore the stan- 
dard mixing length theory was used to evaluate the degree 
of superadiabaticity. In the present paper outer boundary 



conditions will be evaluated by adopting suitable non-grey 
model atmospheres. 

To our knowledge, the most updated model atmosphe- 
res presently available are the ones computed by Brett 
(1995 a,b, hereinafter B95) and the "next generation" of 
the Allard & Hauschildt models (AUard & Hauschildt 1997, 
hereinafter NG97). Both sets of models include updated 
(but not identical) line lists and in both sets the numer- 
ous atomic and molecular opacity sources have been mod- 
eled with the accurate opacity sampling technique (for a 
detailed discussion on this point see also Bessell 1995). At 
the time when present computations have been performed 
NG97 model atmospheres were under-computing (Leggett 
et al. 1996, Chabrier et al. 1996) and results were available 
for solar metallicity only. B95 models are actually avail- 
able from solar metallicity down to metallicities as low as 
[M/H] ~ —2.0. In a recent paper Chabrier et al. (1996, see 
also Baraffe & Chabrier 1996) have discussed the compari- 
son between stellar models computed alternatively adopting 
boundary conditions from NG97, from B95, or from "old" 
(Base) model atmospheres by Allard & Hauschildt (1995, 
hereinafter AH95), emphasizing as a final result "the ex- 
cellent agreement between the observation and the stellar 
models based on B95 and NG97, while models based on 
AH95 clearly underestimate the flux in the V band". As 
pointed out by these authors, this occurrence has to be re- 
garded as a plain evidence of the accuracy recently achieved 
in computing model atmospheres for M dwarfs. 

In order to investigate VLM stellar structures in this pa- 
per we will adopt B95 atmospheric models, which cover the 
ranges AOOOK >Teff > 2mQK, with log g = 4.5, 5.0 and for 
metallicities Z = 0.0002, 0.002, 0.02, 0.04, adopting for the 
mixing length I = 1.5Hp, where Hp is the pressure scale 
height. However, low metallicity VLM models reach larger 
effective temperature (see, e.g.. Paper I) and the "critical" 
temperature (T^Jj?) defining the lower limit for the validity 
of the Eddington approximation also increases above iOOOK 
(see Baraffe et al. 1995 and reference therein). Therefore, at 
the lower metallicities B95 models will be implemented with 
Kurucz's (1993) model atmospheres which are available for 
effective temperature Tg// ^ 3500-fS' and in a large range 
of gravity. One has to notice that Kurucz models lack the 
important H2O opacity contribution which can be of ma- 
jor relevance in cool VLM objects. However, Brett (1995a; 
but see also Bessel 1995) has already found that, regard- 
less of the missing H2O opacity and the different mixing 
length parameter (Kurucz adopts I = 1.25Hp), at 4000i<r 
the two model atmospheres are almost identical, whereas 
at Tgff = 3500if significant differences appear. Thus Ku- 
rucz models can be safely used only for temperatures above 

When using model atmospheres as boundary conditions 
one has to remind that in stellar interiors the radiative flux 
is estimated by means of the diffusion equation which is the 
limit of the transfer equation for large optical depth (see 
Mihalas 1978 for a detailed discussion on this topic). One 
should thus pay attention that the assumed limit between 
the internal structure (where the diffusive approximation 
is used) and the outer atmosphere (where this approxima- 
tion is no longer valid) is located at a Rosseland optical 
depth Tpos large enough for the diffusion approximation to 
be valid. Morel et al. (1994) have recently shown that for 
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Figure 1. The location in tlie HR diagram of solar mctallicity, 
VLM stellar models computed under various assumptions (as la- 
beled) for the boundary condition. Dotted lines connect stellar 
models with the same mass. 

solar abundances a suitable value for Tn^g is about 10. In 
the ATLAS9 code (Kurucz 1993), the monochromatic flux is 
computed through the diffusion approximation when tu (i.e. 
the optical depth at the frequency v) is larger or equal to 20. 
Thus a safe procedure is to connect the ATLAS9 model at- 
mosphere with a model of stellar interior at an optical depth 
insuring the condition Tu > 20 for all frequencies, that it 
means r^os > 20. According to the evidence that physi- 
cal mechanisms related to the energy transport in the outer 
layer are accounted for in a more realistic treatment in the 
atmosphere codes and also to fulfill the quoted condition, 
we fixed the bottom of the atmosphere at thos = 100- The 
same value for the boundary limit has been adopted also by 
Baraffe & Chabrier (1996) and this occurrence should allow 
a safe comparison between our models and the ones provided 
by the quoted authors. 

Nevertheless, it appears interesting to investigate how 
far the evolutionary behavior and in particular the location 
of the models in the HR diagram could be affect by different 
choices concerning the outer boundary limit. For this aim 
stellar models with solar metallicity have been computed 
under boundary conditions taJsen at different values of the 
Rosseland optical depth. 

The computed stellar models cover the mass range 
0.095Mq < A/ < 0.6Mq, assuming an original Helium 
abundance as given by Y=0.27. To be consistent with evo- 
lutionary computations shown in Paper I, we have adopted 
for the mixing length I = 2.2Hp. The adopted value is 
different than the values adopted by Brett (1995) and Ku- 
rucz (1993) in computing their model atmospheres. How- 
ever, it has been verified performing additional numerical 
tests that reasonable variations on the value of / have no 
effects at all on the determination of the stellar radius of 
stars less massive than about 0.5Mq (sec Paper I and refer- 
ences therein), and only minor changes for stars with mass 
around O.GMq. Since in present work, we are mainly dealing 



Figure 2. The HR diagram location of lOGyr old models for 
the chemical composition: Z=0.002 - Y=0.23, computed under 
different assumptions about the outer boundary conditions. 

with stellar models with mass M < 0.6Mq, the results are 
not affected by the choice on the mixing length parameter. 
Moreover, as far as concerns the model atmosphere compu- 
tations, Brett (1995) has clearly shown that a variation of 
the mixing length parameter in the range l.O to 2.0 produces 
only minor effects at effective temperatures around 40007^, 
and that these structural effects decreases rapidly decreasing 
the effective temperature of the models. 

Figure 1 show the HR diagram location of lOGyr old 
stellar models computed adopting r/vo.?: 0.1, 1.0, 100. In all 
models the outer boundary conditions have been obtained 
from B95 model atmospheres alone. For the aim of com- 
parison, the T(r) models as given in Paper I axe also pre- 
sented. As expected, T{t) models axe systematically hotter 
and brighter than B95 models with the same mass. The 
maximum departure occurs at logTgj^ ~ 3.52 which corre- 
sponds to M « O.3OM0, with differences ATe// « 150K 
and AlogL/L0 rj 0.08. The origin of such behavior is re- 
lated to both the H2 dissociation mechanism and to the 
resulting penetration of the convection in the optical thin 
outer atmosphere (Auman 1969; Dorman et al. 1989; Bur- 
rows et al. 1993). 

Again as expected, one finds that at smaller values of r 
the location of the models is largely sensitive to the adopted 
value of the Rosseland optical depth where outer boundary 
conditions are fixed. In fact by changing r/jos &oni 0.1 to 
1.0 the models become cooler by « (70— 80)i^ and fainter by 
AlogL/L0 w 0.05 — 0.1. However, for Tfios > 1-0 the mod- 
els are not significantly affected by changes in the adopted 
value of Tjiog. As a result, one finds the evidence that the 
evolutionary results are not dependent on the choice about 
the point where the atmosphere is matched with the enve- 
lope, once that a sufficiently large value (r^os > 20 — 30) is 
adopted in the evolutionary codes. 

Bearing in mind such a scenario, stellar models for lower 
metallicities have been computed under the above quoted 
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Figure 3. As in Figure 2, but for Z=0.0002 - Y=0.23. 

assumptions, for Z=0.002 and 0.0002 (as in B95), assum- 
ing now Y=0.23 and adopting outer boundary conditions 

either from B95 or from Kurucz's model. For each given 
metalhcity, a set of stellar models still relying on the Ed- 
dington approximation has been also computed in order to 
investigate the effect of incorporating model atmospheres 
in stellar computations as well as to obtain indication on 
the value of the critical temperature T^jf- Figures 2 and 
3 show the results of tfiesc computations. As expected, 
even a quick inspection of those figures shows that at effec- 
tive temperatures around 4000K a good match is achieved 
between the different models, supporting the extension of 
B95 with Kurucz models. More in details, when Z=0.002 
one may notice in Figure 2 the fine smooth transition be- 
tween the T(t) stellar models and the K93 ones which oc- 
curs at Teff w 4AQQK for a w O.SSM© model. The same 
smooth transition occurs between K93 and B95 models at 
T^jf ^ 4000/s: for a ~ O.42A/0 object. This finding ap- 
pears in fine agreement with the results by Brett (1995a, his 
figure 8) when comparing his model atmospheres with the 
ATLAS9 results. The evidence that for T^fj < mOOK the 
K93 model sequence "converges" towards the T{t) models 
is simply due to the missing H2O opacity in the Kurucz's 
(1993) model atmospheres. The stellar models computed 
by adopting the Eddington approximation are hotter and 
brighter than B95 models. The larger discrepancy is present 
at Tgff ~ 38007^ where it is ~ 140A' in effective tempera- 
ture and AlogL/L0 « 0.055 in luminosity. 

The HR diagram location of the stellar models with 
metallicity Z=0.0002 is displayed in Figure 3. In this case 
the match between the T(r) and K93 models is achieved at 
T,.y j ~ 5000-ftr. This behavior appears in good agreement 
with result by Burrows et al (1993) concerning the increase 
of Tgy j when the heavy elements abundance is decreased. 
As expected, even at this very low metallicity B95 models 
agree with K93 ones at T^ff — iOOOK. The maximum dis- 
crepancy between T{t) models and stellar structures com- 
puted by adopting accurate model atmospheres appears at 



Te// « (4100 — 4200)7^, where one finds difi'erences of the 
order of « 180K in T^ff and w 0.003 in L/L© 

3 STELLAR MODELS. 

For each given metallicity, we selected as the "best" 
sequence of models the one obtained by adopting in the 
proper range of validity the T(r), the K93 and the B95 stel- 
lar models, paying attention that a fine and smooth match 
is obtained between the models computed under difi'erent 
assumptions about the outer boundary conditions. Figure 4 
displays the run of present "best" MS in the HR diagram for 
the labeled assumptions about metallicity. Tables 1 through 
3 give the luminosity, the effective temperature, the abso- 
lute visual magnitude and predicted colors (in the standard 
Johnson- Cousins system for VRI and the CIT system for 
K) for lOGyr old "best" models for the various selected 
metallicities. Magnitudes and colors have been evaluated 
by adopting bolometric correction and color temperature 
relation from Kurucz (1993) implemented at efi^ective tem- 
peratures lower than 40007^" with similar evaluations given 
by AUard & Hauschildt (1995) or, for solar metallicity, by 
AUard et al. 1997 (AUard 1996), both properly shifted to 
overlap at each metallicity Kurucz's evaluation at the fitting 
point Te// = 4000if 

The evolution of theoretical predictions concerning the 
HR diagram location of VLM stellar models and the mass- 
luminosity relation has been already discussed by Baraffe et 
al. (1995, hereinafter BCAH95), Baraffe & Chabrier (1996, 
hereinafter BC96) and Chabrier et al. (1996, hereinafter 
CBP96). BCAH95 presented a comparison between stellar 
models based either on the T(r) or on the "old" genera- 
tion of model atmospheres by AUard & Hauschildt (1995). 
CBP96 discussed the effect of different treatments of the at- 
mosphere on the mass luminosity relation, presenting evolu- 
tionary computations which for the lower metallicities still 
rely on the AH95 model atmospheres. When revising this 
paper, improved models for metal-poor low-mass stars have 
been presented by Baraffe et al. (1997, hereinafter BCAH97) 
and for solar metallicity by Chabrier & Baraffe (1997, here- 
inafter CB97), so we have performed some additional com- 
parisons with those updated models. 

Figure 5 compares the HR diagram location of solar 
metallicity VLM models computed by relying either on B95 
or on the Eddington approximation with models presented 
by BCAH95,CBP96 and CB97. One finds that present re- 
sults distribute in between the "old" (BCAH95) and the 
"new" generation models (CBP96) presented by the group 
of Baraffe and coworkers. This is a rather surprising result, 
since present models and CB97 should have quite a simi- 
lar input physics (following the results discussed by BC96), 
and it is not clear to us where the difference is coming from. 
Nor the origin of the (small) differences between BC96 and 
CB97 has been till now discussed in the literature. One also 
finds the curious evidence that most updated CB97 models 
appear in close agreement with our models computed with 
T(r) boundary condition. 

The main difference between the NG97 models (AUard 
& Hauschildt 1997, Lcggctt et al. 1996, CBP96) and the 
B95 ones is the adoption of different lino lists for TiO and 
H2O (by Jorgensen 1994 in NG97 and by Plez et al. 1992 in 
B95). However, BC96 (and also CBP96) compared stellar 
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Figure 4. The HR diagram location of lOGyr old models for the 
labeled assumptions on the chemical composition. 



Table 1. Mass, luminosity, effective temperature, absolute visual 
magnitude and colors for stellar models with solar metallicity and 
Y=0.27, at age equal to lOGyr. 



M/Mq 


log L/Lq 


logTe 


My 


(V-I) 


{V-R) 


{V-K) 


.600 


-1.152 


3.590 


8.954 


1.965 


1.057 


3.817 


.550 


-1.313 


3.570 


9.656 


2.192 


1.147 


4.180 


.500 


-1.468 


3.554 


10.276 


2.352 


1.220 


4.450 


.450 


-1.616 


3.541 


10.806 


2.455 


1.273 


4.635 


.400 


-1.751 


3.531 


11.263 


2.531 


1.311 


4.771 


.350 


-1.881 


3.524 


11.674 


2.587 


1.339 


4.866 


.300 


-2.002 


3.518 


12.056 


2.639 


1.365 


4.953 


.280 


-2.061 


3.516 


12.230 


2.657 


1.373 


4.983 


.250 


-2.159 


3.511 


12.553 


2.707 


1.397 


5.066 


.200 


-2.360 


3.501 


13.254 


2.841 


1.457 


5.285 


.180 


-2.459 


3.495 


13.652 


2.944 


1.502 


5.452 


.150 


-2.638 


3.482 


14.420 


3.152 


1.600 


5.809 


.120 


-2.887 


3.459 


15.721 


3.571 


1.813 


6.557 


.100 


-3.147 


3.429 


17.462 


4.213 


2.175 


7.737 


.098 


-3.181 


3.425 


17.704 


4.301 


2.229 


7.905 


.097 


-3.200 


3.422 


17.868 


4.365 


2.268 


8.029 


.095 


-3.239 


3.419 


18.084 


4.431 


2.309 


8.155 



models computed with their evolutionary code using alter- 
natively the B95 model atmospheres and the NG97 ones, 
obtaining only negligible differences. As discussed by Bessel 
(1995, 1996), one has to remind that for temperatures cooler 
than 3000 A" B95 models have much too strong H2O bands, 
indicating that the H2O opacity is overestimated. Such oc- 
currence should have the effect of decreasing the B95 fits 
of spectra of the coolest M dwarfs by about WOK (Bessel 
1995), the consequences on the HR diagram location of the 
models being correctly valuable only when Brett model at- 
mospheres with updated H2O opacities will bo available. 
However, the differences between present and CB97 models 
largely occur at effective temperatures larger than SOOOif , 
thus the discrepancies can bo hardly ascribed to such an ef- 
fect. Here we can only conclude that such differences deserve 
further investigations, data in Figure 5 giving an indication 



Table 2. As in Table 1, but for metaUicity Z=0.002 and Y=0.23. 



M/Mq 


log L/Lq 


logTe 


Mv 


{V-I) 


{V-R) 




-K) 


.700 


-.514 


3.717 


6.273 


.952 


.540 




1.929 


.600 


-.833 


3.676 


7.240 


1.155 


.661 




2.381 


.550 


-1.030 


3.643 


7.945 


1.379 


.800 




2.788 


.500 


-1.226 


3.617 


8.683 


1.617 


.923 




3.158 


.450 


-1.407 


3.600 


9.292 


1.764 


.997 




3.381 


.420 


-1.507 


3.592 


9.571 


1.797 


1.013 




3.431 


.400 


-1.568 


3.588 


9.739 


1.814 


1.021 




3.458 


.380 


-1.624 


3.584 


9.901 


1.835 


1.029 




3.490 


.350 


-1.704 


3.580 


10.121 


1.855 


1.037 




3.522 


.320 


-1.773 


3.577 


10.312 


1.872 


1.043 




3.549 


.300 


-1.826 


3.574 


10.466 


1.891 


1.048 




3.579 


.250 


-1.981 


3.567 


10.910 


1.941 


1.060 




3.657 


.200 


-2.192 


3.556 


11.553 


2.037 


1.075 




3.806 


.180 


-2.297 


3.549 


11.888 


2.097 


1.085 




3.900 


.150 


-2.486 


3.536 


12.522 


2.229 


1.111 




4.101 


.120 


-2.752 


3.511 


13.546 


2.516 


1.177 




4.533 


.110 


-2.879 


o.4ycs 


14.098 


2.699 


1.230 




4. /y / 


.100 


-3.066 


3.474 


15.087 


3.092 


1.373 




5.370 


.099 


-3.092 


3.170 


15.262 


3.167 


1.405 




5.499 


.098 


-3.119 


3.466 


15.439 


3.240 


1.437 




5.626 


.096 


-3.182 


3.456 


15.865 


3.429 


1.525 




5.921 


.095 


-3.216 


3.451 


16.081 


3.523 


1.572 




6.059 


.093 


-3.294 


3.438 


16.651 


3.793 


1.710 




6.450 


.092 


-3.338 


3.430 


16.991 


3.959 


1.799 




6.683 


Table 3. As in 


laDie 1. 


, but for metallicity Z=0.0002 


and 


Y=0.23. 
















M/Mq 


log L/Lq 


logTe 


My 


{V-I) 


{V-R) 




-K) 


.700 


-.385 


3.758 


5.910 


.805 


.451 




1.580 


.600 


-.743 


3.712 


6.883 


.988 


.553 




1.966 


.500 
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of the degree of freedom still existing in theoretical predic- 
tions. 

As far as the mass-luminosity relation is concerned, 
CBP96 have already discussed the effect of non-gray model 

atmospheres on the reliability of m-L relation, and that dis- 
cussion will not be repeated here. Moreover, Kroupa & Tout 
(1997) (but see also von Hippel et al. 1997) have recently 
presented - still relying on T(t) stellar models - an investi- 
gation on the metallicity dependence of the theoretical mass 
- magnitude relation. Figure 6 shows the m-My and m-Mj 
relations for solar metallicity VLM objects, as obtained by 
using T(t) or B95 atmospheres. For the sake of compar- 
ison we report in the same figure similar predictions from 
BC96, with magnitudes derived from the published lumi- 
nosities according to the above quoted procedure. Figure 7 
shows the mass-luminosity relations in selected photometric 
bands for the "best" models at the two lower metallicities 
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Figure 5. The location in the HR diagram of present stellar 
models for solar metallicity compared with similar models but 
from BCAH95, BC96 and CB97. 



Figure 6. m - My and m - Mj relations for solar metallicity, 
lOGyr old models, for various assumptions concerning the outer 
boundary conditions (see text). 



investigated in the present work. As a relevant point, one 
finds that the m-Mx relation appears scarcely affected by 
metallicity effect, at least for M > (0.11 - O.12)M0, i.e. for 
Mfc < Q.Omag. Such occurrence could be of some help when 
planning observational surveys devoted to investigated the 
mciss function for field stars, whose metallicity is usually not 
well known. 

In panel c) of the previous Figure, we report also 
the semi-empirical m-Mj relation given by Fahlman et al. 
(1989) for a metallicity Z « 0.0001. One finds that such an 
estimate appears in rather good agreement with present the- 
oretical predictions for m-Mj at Z=0.0002 for M > O.IGMq 
i.e. Mj < lOmag. At larger magnitudes, the semi-empirical 
relation crosses our predictions for Z=0.0002, predicting a 
larger magnitude for a given stellar mass. However, at 
this faint end of the m-Mj relation, Fahlman et al. (1989) 
adopted theoretical models by D'Antona (1987), computed 
by using "old" physics both for the EOS and the low temper- 
ature opacity, possibly affecting the semi-empirical results 
for Ml magnitudes above « (9 — 10)mag. 

Since the first derivative of the mass-magnitude relation 
is a key tool in interpreting observed luminosity functions 
in term of the mass function, previous results for the vari- 
ous bands have been best fitted to obtain simple analytical 
relations of the form: 



logM/M© =ao + ai- Mx+a2- M^ + as- M^ + a4-M^ 

Table 4 shows the values of the coefficients for the different 
metallicities together with the value of the standard devia- 
tion a for each relation. 

4 THEORY VERSUS OBSERVATIONS. 

In Paper I we found that a rather satisfactory agreement 
between observational data and theoretical VLM models can 



Figure 7. As in Figure 6 but for Z=0.002 and Z=0.0002. The 
mass- luminosity relation is shown only for "best" models in se- 
lected photometric bands. Panel c) shows also the m-Mj relation 
from Fahlman et al. (1989). 

be achieved even by relying on models based on the Ed- 
dington approximation. Let us here compare present "best" 
models with observations. 

The most relevant observational sample is obviously 
represented by recent Hubble Space Telescope data for lower 
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Table 4. Coefficients of the polynomial regression for the mass - magnitude relation, for the various metallicities and photometric bands 
adopted in the present work. The last column lists the standard deviation. 

Z Mx ctQ ^1 ^2 ^3 ^4 ^ 
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Figure 8. (V, V-I) CM diagram for the lower main sequence 
of NGC639T (Cool et al. 1996) as compared with theoretical 

isochroncs for [M/H]=-2.04 , and for the ages 10, 12 and 14Gyr 
(Cassisi ct al. 1997) shifted to account for a cluster distance 
modulus and reddening (m - M)y = 12.50 and E(V - I) = 0.20. 
The MS locus for VLM structures, for [M/H]=-1.04 t=10 Gyr is 
also shown. 

main sequences in galactic globular clusters. Several CM di- 
agrams have boon already presented (NGC6397: Parcscc, Do 
Marchi & Romaniello 1995, Cool, Piotto & King 1996 and, 
Mould et al. 1996; 47Tuc, M30: King, Cool & Piotto 1996, 



Piotto, Cool & King 1997; NGC6752: Ferraro et al. 1997; 
NGC6656: De Marchi & Paresce 1997). The most tight se- 
quence of VLM stars in a GC appears the one presented 
by Cool et al. (1996), which represents a fundamental tool 
to test the theory of VLM structures. In Paper I it has 
been shown that a largely satisfactory agreement has been 
achieved between our T(t) models and observation. In fig- 
ure 8, we perform the same comparison but using present 
best VLM models for metallicity [M/H]=-2.04 and -1.04, 
adopting from Alcaino et al. (1987) a cluster distance mod- 
ulus (m - M)y = 12.50 and a reddening EiV - I) = 0.20. 
Data in figure 8 have been implemented at larger luminosi- 
ties with isochroncs for the ages 10, 12 and 14Gyr (Cassisi 
et al. 1997), as computed by adopting the same opacity 
evaluations used in this work but the OPAL equation of 
state (Rogers, Swenson & Iglesias 1996) to allow the re- 
quired match with the VLM sequence (see Brocato, Cassisi 
& Castellani 1997a for a discussion on that matter). 

One finds that observational data agree fairly well with 
present theoretical predictions for metal poor models. In- 
teresting enough, one can notice that not reasonable fitting 
can be achieved with the moderately metal rich sequence 
shown in the same figure. Thus theoretical results appear 
in good agreement with current estimates for the cluster 
metallicity, namely [M/H]~ —1.61, where an enhancement 
of a elements by [a/Fe] = 0.30 has been taken into account 
(see Brocato, Cassisi & Castellani 1997b for more details). 
As shown in the same Figure, some residual discrepancies 
between theory and observations still exist, to be eventually 
better understood but only when updated theoretical color 
- Tf,ff relations suitable for metal poor stars will become 
available (AUard et al. 1997). 
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Figure 9 shows the most complete presently available 
CM diagram for stars with known parallaxes, as obtained 
implementing the sample provided by Monet (1992) with 
recent data by Dahn et al. (1995). In Paper I, it has been 
already shown that metal poor T(r) models and BCAH95 
models for [M/H]=-1.5 appear in rather good agreement 
with the location of the hotter subdwarfs sequence. How- 
ever, that paper also disclosed that all theoretical models, 
including BCAH95, failed in accurately reproducing the lo- 
cation of fainter objects for the cooler sequence of stars, 
usually interpreted as the sequence of VLM stars with solar 
metallicity. To test if the use of more accurate outer bound- 
ary conditions can help in reducing the discrepancy between 
observations and theory, the same Figure 9 gives the loca- 
tion of "best" models, for the three metallicities adopted in 
this work, together with the T(r) sequence of models for 
solar metallicity from Paper I and the BCAH95 and BC96 
models (all for solar metallicity). 

An inspection of the figure leads to the following con- 
clusions: 



i) "best" metal poor models rank very well along the hot- 
ter subdwarfs sequence, supporting the indication given in 

Paper I about the lower limit for the metallicity of disk 
subdwarfs and the evidence that the CM diagram location 
of VLM stars appears as a metallicity indicator of unusual 
sensitivity (Paper I, Brocato, Cassisi & Castellani 1997a, 
1997b); 

ii) B95 models significantly improve the fit of the metal rich 
sequence in comparison with oJd T{t) models. Present mod- 
els appear also in best agreement with observation with re- 
spect to BC96 models in the color range 2.2 < {V — I) < 3.0, 
where the main sequence location is strongly affected by 
the adopted treatment for the outer boundary conditions 
(as discussed in Paper I and in BCAH95). Nevertheless, 
one finds that a significant discrepancy still exists for colors 
{V - I) > 3.0 mag, i.e. My > 14 mag; 

iii) Curiously enough, BCAH95 models for solar metallicity 
seem to match the location of the cooler sequence in the CM 
diagram better than both present "best" and BC96 models. 
Since both present and BC96 models have been computed by 
adopting a treatment of the atmosphere much more accurate 
than in BCAH95, this occurrence has to be perhaps regarded 
as an evidence that some other "ingredient" (opacity?, color- 
temperature relation?), adopted in the computations, needs 
further improvements. 

Let us finally compare our VLM stellar models with ob- 
servational data recently provided by Legget et al. (1996), 
who investigated 16 red dwarfs providing bolometric lumi- 
nosities and effective temperatures. These authors already 
found a good agreement between their data and BCAH95 
models. Figure 10 now compares Legget et al. (1996) data 
with present "best" and BCAH97 predictions for various as- 
sumptions about the metallicity, and with BCAH95, BC96 
and CB97 solar metallicity models. Inspection of the Figure 
shows that within the accuracy of observational data one 
may only conclude for a general agreement between theories 
and observation, without making a choice among the dif- 
ferent approaches. In this context, it is worth noting that 
the object Gl 65AB which tends to be cooler than the solar 
metallicity sequence is an unresolved binary system. If this 



Figure 9. My versus {V — I) diagram for faint stars with 
known parallaxes as provided by Monet et al (1992) or Dahn 
et al. (1995). Theoretical predictions from this paper for the 
[M/H]=0.0, -1.04 and -2.04, are also displayed together with the 
theoretical models by BCAH95, BC96 and the T(r) models for 
solar metallicity from Paper I. 



Figure 10. The HR diagram location of VLM object from 
Leggett et al. (1996) as compared with theoretical predictions 
from the present paper or BCAH97 for the various metallicities, 
as labeled. Solar metallicity models from BCAH95 and BC96 are 
also displayed. 

system would resolved into two similar components, these 
would lie near the point marked by an arrow (Leggett et al. 
1996), in fine agreement with the theoretical prediction. 

Figure 11 compares theoretical m-My relations from 
the present paper and from BCAH95 and BC96og0 with ob- 
servational data for nearby binaries by Henry & McCarthy 
(1993). The two objects at My ~ 8.6 mag (fisted as Gl 
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Figure 11. Mass-luminosity relations for the labeled assump- 
tions on the stellar metallicity, as obtained in this work, by 
BCAH95 and BC96. Observational data have been provided by 
Henry & McCarthy (1993). 



Figure 12. Mass - absolute visual magnitude relations for solar 
metallicity stellar models. Present results are plotted as solid 
line, models by CB97 as a dashed curve. Models of BCAH95 and 
BC96 are plotted as open squares and open triangles, respectively. 
Observational data correspond to the empirical m — My relation 
derived by Kroupa, Tout & Gilmore (1993). 



677A and Gl 677B) which deviate significantly from theo- 
retical prescriptions, belong to the same binary system, with 
no large accuracy of the orbital parameters since the orbit 
has been followed for less than a revolution (Henry & Mc- 
Carthy 1993) Again one finds that theory and observations 
appear in satisfactory agreement, without allowing a choice 
among the diflerent theoretical approaches. 

Figure 12 finally compares the mass luminosity function 
given by Kroupa ct al (1993) with present and previous com- 
putations on the matter (for solar metallicity), as labeled. 
Now one finds that the most recent results appear in better 
agreement, without any clear indication allowing a choice 
between present models and CB97. In the above quoted pa- 
per (and reference therein) Kroupa ct al. discussed also the 
relevance of the run with My of the derivative of the mgiss- 
luminosity function dm/dMy. Therefore, we have decided 
to compare the derivative of the mass - luminosity function 
obtained by adopting the most updated VLM stellar mod- 
els, presently available, with the derivative of Kroupa et al.'s 
(1993) relation, as given by Kroupa & Tout (1997). Figure 
13 shows that present models foresee a maximum in the 
derivative whose location appears in better agreement with 
observational prescriptions than CB97 models do. Thus in 
this respect it appears that present models work better. 

5 CONCLUSIONS. 

In Paper I, we investigated the effects of new physical in- 
puts as equation of state and low temperature opacities on 
the evolutionary properties of VLM stellar objects, by com- 
puting stellar models relying on the Eddington approxima- 
tion for the treatment of outer boundary conditions. In this 
work, we devoted our attention to the effect of improving the 
treatment of the atmosphere, by adopting Brett's (1995a,b) 
non-grey model atmospheres. 



Figure 13. The derivative of the m — My relation as function of 
My. The solid curve is the derivative of present models and the 
dashed line corresponds to the derivative obtained adopting the 

CB97 stellar models. The observational data correspond to the 
derivative of the mass - luminosity relation provided by Kroupa 
& Tout (1997). 

Present models have been compared with other VLM 
stellar models in the current literature, as given by BCAH95, 
BC96, BCAH97 and CB97. Comparison with observational 
data for metal poor GC stars or solar neighborhood dwarfs 
shows a reasonable agreement. This result can be regarded 
as a plain evidence that the solution of the long-standing 
discrepancy between theoretical predictions and observed 
HR diagram location of VLM stars is no more a tantaliz- 
ing goal. However, a not negligible discrepancy still exist 
between stellar models and CM location of the cooler se- 
quence of VLM objects with known parallEixes, for magni- 
tudes larger than My > 14mag. This occurrence and the 
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presence of still significant uncertainties in the observational 
data - in particular concerning the effective temperature - 
confirms (as pointed out by Legget et al. 1996, Bessell 1995, 

CBP96) that there arc still improvements to bo performed, 
both in the models (structural and atmospheric) and in the 
observations. 
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